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Abstract 

Based on experimental observations in A X M02 (A = Na, Li; M = Co, Ni), a model for suppression 
of magnetic frustration by electron doping in a nearest-neighbour antiferromagnetic triangular 
lattice is presented. It is found that frustration can be quantified, as determined by geometry 
and bond-counting, and its magnitude is a non-monotonic function of x. A mean-field calculation 
provides temperature-dependent magnetization, spin-entropy and heat capacity. Low-doping (x = 
0.25, 0.33) results in a highly frustrated regime. A0.5MO2 has strongest order and no frustration, 
while high doping (x = 0.67, 0.75) leads to low frustration and higher spin-entropy. The results 
agree with experiments including neutron scattering, spin entropy-driven thermoelectricity and ion 
ordering. 

PACS numbers: 81, 75.10.-b, 61.14.-x 
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The discovery of superconductivity in N a^CoC^.y H 2 has reignited the importance of 
electrical and magnetic properties of the doped triangular latticed. The structure con- 
sists of a layered triangular lattice of Co ions within oxygen octahedra, with Na ions and 
water molecules intercalated between these layers^. The Na x Co02 {x ~ 0.6 - 0.8) 
compounds are important materials for thermoelectric applications as they exhibit a large 
thermopower, low resistivity and low thermal conductivity^. Systematic experiments^ ' 6 ' 7 ' 8 ' 9 
and theor y 10 ' 11 ' 12 ' 13 as a function of x have identified charge, spin, and alkali metal-ion (A) 
ordered phases at particular x values in the isostructural family of compounds A X M02 (A = 
Na, Li; M = Co, Ni). These include the material used in present-day rechargeable Lithium 
batteries^. Strong correlations (on-site Coulomb energy ~3-5 eV in the cobaltate a 11 ' 15 ' 16 ) 
are important for spin-charge order, thermodynamics and transport properties. Doping, 
frustration and fluctuations in a triangular lattice (TL) lead to well-known results, includ- 
ing the resonating valence bond state and superconductivit y 17 ' 18 ' 19 . 

Nao.5Co02 undergoes a charge and magnetic ordering metal-insulator transition with 
T^' 20 ' 21 . For x < 0.5 a paramagnetic correlated metal phase exists, with a significantly 
enhanced in-plane magnetic susceptibility y^. In contrast, for x > 0.5 a " Curie- Weiss" (CW) 
metallic state (localized Co spins with a CW-% and a T-linear resistivity) is established^. 
The large thermopower for x ~ 0.6 — 0.8 is due to enhanced electronic spin entropy, that 
gets suppressed upon applying a magnetic field 4 . Interestingly, for x ~ 0.66 — 0.75, x(T) fits 
a CW law with an effective Weiss temperature 9 that increases from ~ -50 K to -125 K— . 
This suggests increasing AF correlations with increasing x, for x > 0.5, and a long-range 
ordered AF state is known for x ~ 0.75, 0.82 22 ' 23 . Surprisingly, neutron inelastic scattering 
indicates ferromagnetic in-plane spin-fluctuations for Nao^CoO?^. A similar result is also 
known for NaNiO^. For Nao.6-o.sCo02, the average local moment per Co site is estimated 
to be larger than that expected in a simple chemical picture of electron doping due to Na 
ions. The doped electrons are expected to form low spin 5 = Co 3+ sites in a background 
of S = 1/2 Co 4+ ions (x = implies all S = 1/2, Co 4+ , while x = 1 implies a lattice of 
S = 0, non-magnetic Co 3+ sites). In particular, the experimentally obtained average local 
moment for x = 0.75 is ~1.1 while expected value is ~0.43 [i^. 

Recent NMR measurements which probe the local spin indicate dynamic spin fluctua- 
tions and interlayer coupling, mediated by the Co-O-Na-O-Co exchange path 26 . Another 
recent NMR study has concluded the existence of 3 types of Co ions in Nao.7Co0 2 27 . More 
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intriguingly, two independent neutron scattering studie d 28 ' 29 on high Na content single crys- 
tals (x = 0.75,0.82), corresponding to a very dilute lattice of S = 1/2 Co 4+ ions, show 
that obtained spin wave dispersions can be accurately modelled only by assuming S = 1/2 
Co 4+ moments at all Co sites. This suggested possible phase separation in high Na-content 
materials. Both these studies obtain a ferromagnetic in-plane (Jn) and an AF out of plane 
(J_i_) exchange interaction of comparable magnitudes, indicative of 3-dimensional magnetic 
fluctuations. The studies highlight several unconventional aspects: (i) Jii is comparable in 
magnitude to Jj_, inspite of a layered 2D structur e) 28 ' 29 (ii) 9 inferred from fitted exchange 
parameters in terms of a local-moment picture, namely -(3J|| + Jj_), is positive, while 9 
obtained from x is negative^. The possibility of intermediate spins (S = 1) for Co 3+ was 
invoked to explain this behavior (iii) The observed sharp spin wave modes along the c-axis 
requires that, in a phase separated scenario, clusters should be aligned along c-axis over 
many layers^ (iv) A Wigner crystal model could be reconciled with observed spin-waves if 
Co 3+ ions have a moment similar in size to Co 4+ ions^. 

Another interesting aspect of A X M02(A = Li, Na; M = Co, Ni) is the alkali metal 
ion (A)-ordering observed in this isostructural family. Electron diffraction (ED) studies 
on Na x Co02 found ordered Na ion-Na vacancy superlattices, beyond the simple hexagonal 
structure^. It identified a novel structural principle: the presence of lines of Na-ions and 
vacancies (along the [110] direction), rather than a simple maximization of Na-Na ions in the 
layered structure. This behavior arises from a minimization of electrostatic energy, combined 
with a constraint of occupying Na ions on discrete lattice sites^ 2 -. There exist 2 types of A 
sites in A X M02: Al (or Nal: above or below the M ions in the 2D triangular layer) and 
A2 (or Na2: above or below centroids of triangles formed by M ions). The Nal site was 
shown to be 67 meV higher in energy than Na2, making Na2 the preferred sitei 2 -. ED studies 
confirm that Na2 sites are preferentially occupied, while Nal sites are dominantly empty, 
for all x between 0.15-0.75. The only exception is the x = 0.5 compound which has Nal 
and Na2 sites equally occupied^. The calculated formation energies as a function of x for 
Li^NiO^ and Na^CoOjr^ show very similar behavior, identifying ion ordering at particular 
x values. 

Given these unconventional magnetic observations and the A-ion ordering in the AJ.M02 
series, we resorted to a simple model which aims to study the relation between doping 
fraction x, A-ion ordering, and temperature dependent magnetization, spin-entropy and 
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heat capacity. We describe the M spins by an antiferromagnetic Ising model on a triangular 
lattice. The doped electrons are also described as Ising spins. While our model applies 
equally to the isostructural A X M02 family, we mainly focus on Na^CoC^ for comparion 
with experiments. 

We make the following ansatz: a doped electron occupies an A2 site and, by symmetry, 
interacts equally with all three Co 4+ sites constituting a host triangle. It does not simply 
transform a single Co 4+ site to a Co 3+ site. This electron is assumed to couple antiferromag- 
netically with the three Co 4+ spins of the host triangle, and this makes the coupling between 
Co 4+ spins of the host triangle to become negligible. Thus, doping results in a local suppres- 
sion of frustration (Figs, la and lb). The ansatz may be viewed as dynamic fluctuations of 
spins due to doping, as concluded by NMR and neutron scattering measurement o 26 i 27 i 28 i 29 . It 
yields a quantification of frustration and consistency with many experiments: neutron scat- 
tering, A-ion ordering, heat-capacity and also provides an explanation for enhanced spin 
entropy in the presence of low frustration. 

We investigate x-values which effectively represent the parent superconducting compo- 
sitions (1/4, 1/3), the charge and spin-ordered compound (1/2), and the thermoelectrics 
(2/3, 3/4). These fractions can be classified on a geometric basis as (i) hexagonal fractions: 
x = 1/3 and 2/3 correspond to 1 and 2 occupied Na2 sites in a hexagonal cell, respectively 
(e.g. Fig. 1(c) for x = 1/3), and (ii) rhombic fractions: x = 1/4, 1/2 and 3/4 correspond to 
occupying 1, 2 and 3 Na2 sites in the rhombic cell, respectively (e.g.Fig. 1(d) for x = 1/4). 
Lowest energy structures for these x- values obtained using the model are shown in Figs. 1(e) 
to l(i). A doped electron from an Na2 site aligns the three Co spins of the host triangle 
antiparallel to it. The remaining Co spins are also antiferromagnetically coupled Ising spins, 
aligned so as to minimize the total energy of the unit cells. The deduced Na-ion ordering for 
the lowest energy unit cells (red dotted lines in Figs. l(e)-l(i)) match the experimentally ob- 
tained ion ordering for all x, except for x = 1/2. More significantly, it yields a quantification 
of frustration, defined as the number of unsatisfied bonds (blue dashed lines in Figs. 1(e) 
to l(i)). The unsatisfied bonds are bonds aligned ferromagnetically inspite of the antifer- 
romagnetic coupling. The bonds aligned ferromagnetically for host triangles (containing an 
A2 site) do not count as unsatisfied bonds. The fraction f u b{x) is the number of unsatisfied 
bonds normalised by the number of Co spins in a unit cell: f u b{x) = (3/4, 1, 0, 1/3, 1/4) for 
x = (1/4,1/3,1/2,2/3,3/4) (Fig. 3(c)). Since f ub (0) = 1, the fractions 1/4 and 1/3 are 
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FIG. 1: (Color online) (a) Schematic of frustration, (b) ansatz: a doped electron suppresses frus- 
tration by aligning three NN Co spins ferromagnetically, (c) hexagonal cell, (d) rhombic cell, (e)-(i) 
Ordered structures obtained using the ansatz for re = 1/4, 1/3, 1/2, 2/3, 3/4 (red dotted lines: 
lowest energy unit cells, blue dashed lines: unsatisfied bonds). For x = 1/2, black full line shows 
equivalent orthorhombic unit cell 
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FIG. 2: (Color online) Magnetization m^T), for x = 1/4, 1/3, 1/2, 2/3, 3/4 (J' = 1), obtained 
from mean-field theory for distinct sites i in a unit cell. 

considered highly frustrated. The suppression of frustration is found to be a non-monotonic 
function of x. For x = 1/3 (Fig. 1(f)) and 2/3 (Fig 1(h)), we need to use a two hexagon 
unit cell for the lowest energy structures; similarly for x = 3/4, the unit cell consists of two 
rhombii (Fig. l(i)). However, for x — 1/4 (Fig. 1(e)) and 1/2 (Fig. 1(g)), a single-rhombus 
unit cell is appropriate. These unit cells (Figs. l(e)-l(i)) also correspond to the lowest f u b 
values, e.g. a single hexagon unit cell for x = 1/3 (Fig. lc) leads to very high frustration 
with f u b = 2. For x — 1, the analysis shows all Co spins are ferromagnetically aligned with 
/u&(l) = 0; the unit cell in this case is a simple rhombus of two triangles that share a bond, 
with one doped electron per unit cell. It is noted that our model for x = 1 is related to the 
delta-star transformation of Baxter and Wu for the 2D triangular lattice, which was solved 
exactly in terms of three-spin interactions to yield a phase transition 30 . 

The above discussion is at zero temperature, T = 0, with all fractions exhibiting com- 
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plex magnetic long-range order. Since / u & varies non-monotonically with x, the magnetic 
transition temperatures, T c (x) are expected to show non-monotonic behaviour. It is also 
important to investigate T-dependent sub-lattice magnetization, spin-entropy and heat ca- 
pacity. We explore these aspects by performing a finite-temperature calculation using an 
Ising-like Hamiltonian, 



Here, (ij) are nearest neighbour pairs of the unit cell. For proper counting, Jy = J/2 
whenever bond (ij) is shared between two unit cells. Our ansatz implies that coupling 
between Co spins of a host triangle is zero, so that = for all such pairs (ij); and for 
all other pairs (ij), J%j = J- We choose energy units such that ,1 — 1. In the second term 
above, k = 1,2, njv a , n> Na = xrico being the number of Na spins a (nc is the number of 
Co spins s) in the unit cell; i^ runs over the three vertices of host triangle k. The parameter 
J' is the Na-Co coupling. Eqn. 1 thus describes the unit cell Hamiltonian for any doping x. 

We perform a mean-field calculation, making the approximation SiSj ~ rriiSj + rrijSi — 
miirij, where magnetization = (sj). This calculation is straightforward, and we note that 
the doped electron spin, the Co spins neighbouring it, and those not neighbouring it, are 
distinct in mean-field theory. We solve self consistency equations for all the order parameters 
m ; . The number of distinct spins, w(x), depends on the unit cell geometry and doping, and 
determines the high-temperature spin entropy S(x,T > T c (x)) = T\nw(x). We calculate 
T-dependent rrii at each distinct site, as well as spin-entropy S and heat capacity C v for all 
fractions x. These thermodynamic quantities are presented in Figs. 2 and 3. 

Fig. 2 shows magnetization m(T) as a functon of x for J' = 1 which indicates magnetic 
order occurs below a critical temperature T c (x) for all x. The function T c (x) anticorrelates 
with f u b(x), and is also non-monotonic (Fig. 3(c)). The plots also suggest rhombic fractions 
are magnetically more stable than hexagonal fractions: in spite of larger frustration for 
x = 1/4 than for x = 2/3, the T c values are nevertheless comparable. 

It is interesting to see that frustration f u b(%) is maximum for x = 1/3, which shows 
the highest superconducting T sc upon water intercalation^. While we obtain magnetic order 
for x = 1/3, it has the lowest magnetic ordering T c . For x = 1/4, as well as 2/3, the 
calculated magnetic T' c s are higher and comparable. While there is no report to date of 
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FIG. 3: (Color online) Plots of (a) spin-entropy and (b) heat capacity as function of tempera- 
ture(T/J), for J' = 1). (c) Plots of frustration f u b{x) and magnetic-ordering temperature T c (x) 
showing non-monotonic behaviour. f u b(x) is maximum for x = 1/3. T c {x) anticorrelates with 
/un- 
observed magnetic ordering in Na^OoC^ (x = 1/4,1/3 and 2/3), electrochemical measure- 
ments suggest an order-disorder transition around x = 1/3- . It is also known that Lij-CoC^ 
and Na x CoC>2 show ion-ordering at x — 1/4, 1/3, 2/3 or very close to these values^*- For 
x = 1/3,2/3 and 3/4, the Na-ion and vacancy ordering occurs along the [110] direction 
from our results (Fig 1(f), 1(h) and l(i)). This is in good agreement with ED results for 
x = 0.3,0.64 and 0.73, which show incommensurate ordering 9 , while our results correspond 
to the nearest commensurate ordering. The results also match the ion ordering in x — 1/4 
with a (2a x 26) superlattice^^*^. The obtained ion ordering for x — 1/2 gives a 'Star of 
David' or 'Kagome' structure within the triangular lattic e - 1 ^ 19 . This is not consistent with 
the x = 1/2 observed Na ion-ordeA Recent work shows that for x = 1/2, the formation of 
a divacancy cluster results in occupancy promotion of an Na2 site to an Nal site^. Since we 
have neglected Nal sites in our model, we cannot address the Nal site ion order. Nonethe- 
less, the Co-lattice magnetic structure obtained from our model for x = 1/2 (2a x 2b or 
equivalently, 2a x \/3a ; Fig. 1(g)) matches extended Hubbard model calculations^ and 
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recent neutron diffraction studie d 20 ' 21 , which show chains of antiferromagnetically ordered 
moments. While we obtain alternately anti-ferro and ferro chains parallel to each other, ref. 
21 models it as anti-ferro chains alternating with non-magnetic chains. However, they also 
point out that the S = spin Co 3+ sites which are supposed to constitute non-magnetic 
chains may actually correspond to Co ions with fluctuating moments which do not order, 
consistent with absence of unit-valence charge disproportionatio n 11 ' 21 . The present results 
also show that the magnetic order is strongest for Nao.sCoC^, as is known experimentally, 
and frustration is completely suppressed ( f u b(x = 1/2) = 0). For the larger doping fractions, 
x = 2/3, 3/4, the frustration is less than for x = 1/3, 1/4, but is still not zero. 

Fig. 3 presents (a) spin-entropy(S') and (b) heat capacity(C„) as functions of tempera- 
ture^/ J). The plots show that C v for x = 2/3 and 3/4 are anomalous (region marked by 
arrow in Fig. 3b). This anomalous behavior has been experimentally observed in the elec- 
tronic part of C v for x = 3/422. Also, Fig. 3a shows that changes in S are higher for x = 2/3 
and 3/4 compared with x = 1/4, 1/3 and 1/2. This is the doping regime which experi- 
mentally shows enhanced thermoelectricity that can by suppressed by an applied magnetic 
field, indicating a magnetic origin. Our results thus confirm the same. We investigated this 
anomaly for x = 3/4 as a function of J', and found that this behaviour shows up only for 
J' > 0.7. If we associate J' and J to be m Jj_ and Jy obtained from neutron scatterin g 2 ^' 2 ? , 
our results are consistent with the experimentally estimated J± ~ (0.73 — 2.1) Jn. 

In summary, we have presented a model for magnetism and related ion ordering for 
the A X M02 (A = Li, Na; M = Co,Ni) family of compounds, in which frustration and 
geometry play an important role. A mean-field calculation provides doping and temperature- 
dependent magnetization, spin-entropy and heat capacity. Low doping (x = 0.25, 0.33) 
results in a highly frustrated regime. A0.5MO2 has strongest order and no frustration, while 
high doping (x = 0.66,0.75) leads to low frustration and higher spin-entropy in the model. 
The obtained results are in general agreement with neutron scattering, spin entropy-driven 
thermoelectricity and ion-ordering. 
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